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Abstract 

(February 1, 2008) 

We represent QCD at the hadronic scale by means of an effective Hamiltonian, 
i7, formulated in the Coulomb gauge. As in the Nambu-Jona-Lasinio model, 

o . 

. chiral symmetry is explicity broken, however our approach is renormalizable 

CN 

^ . and also includes confinement through a linear potential with slope specified 

c3 



by lattice gauge theory. This interaction generates an infrared integrable sin- 
gularity and we detail the computationally intensive procedure necessary for 
OO . numerical solution. We focus upon applications for the u, d, s and c quark 

fiavors and compute the mass spectrum for the pseudoscalar, scalar and vec- 
tor mesons. We also perform a comparative study of alternative many-body 

o , 

. techniques for approximately diagonalizing H: BCS for the vacuum ground 

JL . state; TDA and RPA for the excited hadron states. The Dirac structure of 

D , 

r~| , the field theoretical Hamiltonian naturally generates spin-dependent interac- 



X 



tions, including tensor, spin-orbit and hyperfine, and we clarify the degree of 
level splitting due to both spin and chiral symmetry effects. Significantly, we 
find that roughly two-thirds of the n-p mass diff'erence is due to chiral sym- 
metry and that only the RPA preserves chiral symmetry. We also document 
how hadronic mass scales are generated by chiral symmetry breaking in the 
model vacuum. In addition to the vacuum condensates, we compute meson 
decay constants and detail the Nambu-Goldstone realization of chiral sym- 
metry by numerically verifying the Gell-Mann-Oaks-Renner relation. Finally, 
by including D waves in our charmonium calculation we have resolved the 
anomalous overpopulation of J/\l' states relative to observation. 

PAC number(s): 12.39.Pn, 12.40.Yx 



I. INTRODUCTION 



In a series of publications |T|-^ an ambitious QCD program has been initiated to compre- 
hensively investigate hadron structure. The theoretical formulation entails renormalization 
and utilizes established many-body techniques to approximately diagonalize an effective con- 
fining Hamiltonian. This paper, a detailed continuation of our recent letter , focuses upon 
the quark sector and reports numerical results for mesons complementing our previous gluon 
study [§. 

Over the years there have been many meson investigations, from the early, simple non- 
relativistic constituent quark model calculations to more involved relativistic, field theoreti- 
cal approaches implementing current quarks and spontaneous chiral symmetry breaking. A 
common shortcoming of these analyses is an inability to consistently reproduce the physical 
mass spectrum of the scalar and pseudoscalar mesons. Our paper addresses this issue and 
significantly extends the pioneering work of the Orsay group 0, Adler and Davis 0, and 
the Lisbon investigators [0. In our approach the exact QCD Hamiltonian in the Coulomb 
gauge is modeled by an effective, confining Hamiltonian, H, that is fully relativistic with 
quark field operators and current quark masses. However, before approximately diagonal- 
izing H, a similarity transformation is implemented to a new quasiparticle basis having a 
dressed, but unknown constituent mass. As described in Sec. II, this transformation entails 
a rotation which mixes the bare quark creation and annihilation operators. By then per- 
forming a variational calculation to minimize the ground state (vacuum) energy, a specific 
angle and corresponding quasiparticle mass is selected. In this fashion chiral symmetry is 
dynamically broken and a non-trivial vacuum with quark condensates emerges. This treat- 
ment is precisely analogous to the Bardeen, Cooper, and Schrieffer (BCS) description of 
a superconducting metal as a coherent vacuum state of interacting quasiparticles combin- 
ing to form condensates (Cooper pairs). Excited states (mesons) can then be represented 
as quasiparticle excitations using standard many-body techniques which in this work will 
be the Tamm-Dancoff (TDA) and random phase approximation (RPA) methods. The two 



2 



treatments are truncated at the one quasiparticle, one quasihole level and then numerically 
compared. Our RPA analysis confirms and extends the early work of Ref. which utilized 
an extended Nambu-Jona-Lasinio mean field approach. 

Two other comments are in order before proceeding. First, there are several reasons for 
choosing the Coulomb gauge framework. As discussed by Zwanziger [0, the Hamiltonian 
is renormalizable in this gauge and, equally as important, the Gribov problem (V ■ A = 
does not uniquely specify the gauge) can be resolved (see Refs. for further discussion). 
Related, there are no spurious gluon degrees of freedom since only transverse gluons enter. 
This ensures all Hilbert vectors have positive normalizations which is essential for using 
variational techniques that have been widely successful in atomic, molecular and condensed 
matter physics. Second, due to Fock space truncations our analysis is not Lorentz invariant. 
However, we only plan to use one frame and do not compute hadron form factors with this 
method. Interestingly, violating Lorentz non-invariance implies a prefered reference frame, 
which, as selected by chiral symmetry breaking, is the condensate rest frame. 

This paper consists of six sections and two appendices. The next section introduces our 
effective, QCD inspired Hamiltionian and developes the BCS vacuum treatment leading to 
the quasiparticle mass gap equation. We also compare our approach to the classic Nambu- 
Jona-Lasinio model. In Sec. Ill we detail our numerical, supercomputer solution of the 
gap equation along with the quark condensate and constituent mass values. Sections IV A 
and IV B describe the TDA and RPA, respectively, while Sec. IV C addresses weak decays 
and Sec IV D presents a derivation of the Gell-Mann-Oakes-Renner relation. The TDA 
and RPA meson spectra are compared and discussed in Sec. V. This section also includes 
results from a simple SUf{3) flavor mixing analysis for the 77-77' system and our predictions 
for the charmed mesons. Conclusions and future work are summarized in Sec. VI. Finally, 
Appendix A provides further details regarding the BCS transformation and vacuum state 
while Appendix B presents the most general TDA equation for arbitrary angular momentum. 
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II. HAMILTONIAN AND MASS GAP EQUATION 



A. Effective Hamiltonian 

By introducing a phenomenological confining potential, Vl, the QCD Coulomb gauge 
Hamiltonian for the quark sector can be replaced by an effective Hamiltonian 

H = J ■ V + /3m)^(f) -^J dxdyp''{x)VL{\x - y\)p''{y) (2.1) 

where ^ , m and p°'{x) = '^'^ {x)T°''$ {x) are the current (bare) quark field, mass and color 
density, respectively (for a more complete discussion, especially for the heavy quark sector, 
consult Refs. P,pi|). For notational ease the flavor subscript is omitted (same H for each 
flavor) and the color index runs a = 1...8. Motivated by lattice gauge studies we adopt 
a linear confining interaction, Vl = — y\, with slope a = .18 GeV^ also specified by 
lattice and Regge phenomenology. In our analysis we have also performed calculations 
with and without the leading QCD canonical or Coulomb (one-gluon exchange) interaction, 

2 

= with as = ^ = .4. For most observables, especially the meson mass spectrum, 

the Coulomb interaction is not important and can be omitted. This can be understood by 
noting that in momentum space, where we perform all calculations, the two interactions 
have the same sign, i.e. V{r = \x — y\) = Vc + Vl 



V(k) = / drV{r)e 



-ik-f 

Because the meson wavefunctions have a finite momentum distribution, most static me- 
son properties are predominantly governed by the infrared (IR), or low, momentum region 
where the confining potential dominates. Including the Coulomb interaction is then roughly 
equivalent to using a slightly larger string tension, a. There are certain observables, and in 
particular the gap equation detailed below, for which the Coulomb interaction is ultra-violet 
(UV) divergent. In such cases we regularize with a cut-off parameter and then could renor- 
malize to remove cut-off sensitivity using one of our renormalization procedures detailed in 



Refs. plffl for the gluon sector. In this paper we only present unrenormahzed results since 



this program is still in progress |[TT||T3 and has not yet been completed for the quark sector. 
This is an additional reason for omitting the Coulomb interaction. Hence, with the exception 
of the current quark masses (we use rriu = rrid = 5 MeV, uis = 150 MeV, rric = 1200 MeV), 
our approach entails only one pre-determined parameter which also sets the hadronic scale, 
^ = 424 MeV. 

We further note that even though the confining potential is IR divergent, this singularity 
is cancelled (see Ref. [§) in both the mass gap equation and all calculations for associated 
observables. Hence, the problem is the delicate numerical evaluation of this integrable 
singularity which we discuss in Sec. HI. 

Finally, we stress that in constituent quark models free quarks can exist which requires 
imposing color confinement. However, as demonstrated in Refs. [§-01 the Lorentz structure 
of our Coulomb gauge density-density confining interaction only permits stable solutions for 
color singlet states. Therefore, confinement naturally emerges in our approach. 



B. BCS transformation and gap equation 



We now wish to solve = E"^ as accurately as possible. In this subsection we focus on 
the ground state and introduce the Bogoliubov-Valatin, or BCS, transformation. We begin 
by recalling the plane wave, spinor expansion for the quark field operator 

= E / [u.x{k)Kx{k) + v^,{-k)dlxi-k)] (2.3) 

with free particle, anti-particle spinors McA, Vc\ and bare creation, annihilation operators 
6cA, dcx for current quarks, respectively. Here the spin state (helicity) is denoted by A and 
color index by c = 1, 2, 3 (which is hereafter suppressed). Because we can expand \1/ in terms 
of any complete basis we may equally well use a new quasiparticle basis 



dk 

(M 



Uxik)Bx{k) + Vx{-k)D{{-k) 



Ak-x 



(2.4) 
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entailing quasiparticle spinors Ux, Vx and operators Bx, Dx- The Hamiltonian is equiva- 
lent in either basis and the two are related by a similarity (Bogohubov-Valatin or BCS) 
transformation. The transformation between operators is given by the rotation 



Bx{k ) = cos^bx{k ) - Xsm^d{{-k ) 

Dx{-k ) = cos ^-^dx{-k ) + A sin ^-^b{{k ) 
involving the BCS angle 9^ — 9{k). Similarly the rotated quasiparticle spinors are 



(2.5) 



Ux{k) = cos Y^xik) - A sin ^vxi-k) = 



Vx{-k) = cos^vxi-k) + Xsin^uxik) = 



^Jl + sin (j){k) Xx 



sin (t){k) a-k Xx 



(2.6) 



l-sin0(A;) a ■ k xx 



^1 + sin0(A;) xx 

where xx is the standard two-dimensional Pauli spinor. We have also introduced the gap 
angle, 0^ = which is related to the BCS angle, 9/2, by = 9 + a where a is the 



current, or perturbative, mass angle satisfying sin a = m/Ek with Ek = ywF+W. Hence 



• , n k ■ n 

sm 0fe = — COS 9k + — sm 9k 

COS (f)k = — COS 9k- — sm 9k . 

^k 

Similarly, the perturbative, trivial vacuum, defined by &a|0) = c?a|0) = 0, is related to 
the quasiparticle vacuum, -BaI^) = Dx\^) = 0, by the transformation 



1^^) =exp (-E / ^Atan|6l(^)4(-^)) |0) . 



9, 



(2.7) 



A - (2^ 

In this paper we will denote the BCS vacuum by |Q) (in Sec. IV we introduce the RPA 
vacuum labeled \^Irpa))- Expanding the exponential and noting that the form of the op- 
erator b^d^ is designed to create a current quark/ antiquark pair with the vacuum quantum 
numbers, clearly exhibits the BCS vacuum as a coherent state of quark/antiquark excita- 
tions (Cooper pairs) representing '^^'^^Lj ='^Po condensates. One can regard tan ^ as the 
momentum wavefunction of the pair in the center of momentum system. 
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We now seek an approximate ground state for our effective Hamiltonian by minimizing 
the BCS vacuum expectation, {Q\H\Q). We do this variationally using the gap angle, 0^, 
(not the BCS angle) which leads to the gap equation, 6{Q\H\Q) = 0. After considerable 
mathematical reduction, the nonlinear integral gap equation follows 

2 r dq ^ 

ksm(j)k - mcos(f>k = - J j^^V{\k - ^)[sm(j)kCOS(f>gk ■ q - sincpgcoscpk] ■ (2.8) 

The angular integrals can be analytically evaluated (see Appendix B) to give 

2 1 

k sin (l)k - m cos (l)k = -^^^ g^(ig[cos sin 0fcFi(fc, g) - cosc/)^ sin(/)gVo(fc, ?)] (2.9) 
where 

— IGttct 



(F - g2)2 



k'^q'^ 



k — q J V(^^ ~ I'^y , 

corresponding to the linear potential above. Similar expressions for the Coulomb potential 

are given in Appendix B. 

There are several alternative ways to derive this same gap equation. One is through the 
Ward identites. Another is by requiring cancellation of the anomalous Bogoliubov terms in 
the 2-body part of the newly normal ordered Hamiltonian. The latter is necessary to stabilize 
the vacuum and is also equivalent to minimizing the 0-body constant energy splitting the 
BCS and trivial vacua (see Ref. 0). 

Numerically we actually solve a different form of the gap equation, originally obtained 
by Adler and Davis 0, that is more familiar to the solid state community. They use the 
function ip^ = ip{k) related to our gap angle by 

sm (j)k = ——2 

COS0fc = 2 



with corresponding gap equation 



7 



^^'^-y(i-^.) = ^X YTJl ■ ^'-'^^ 

Examination of Eqs. ( |2.9|j2.1(]| ) reveals that the diverg = g is an integrable 

singularity for the linear potential since for all k the integrands vanish at = g. This is not 
the case for the Coulomb potential UV singularity. However, since it naturally emerges from 
the canonical QCD Hamiltonian (one gluon exchange), we retain the option of including 
this potential for selected calculations and use a cut-off to regulate its ultraviolet (UV) 
divergence. 

The solution of the gap equation (see Sec. Ill) leads to a vacuum quark- ant iquark 
condensate given by 

{qq) = (^]|^(0)^(0)|^]) = -4t / k'^smcPk dk (2.11) 

which is quadratically divergent for non-zero current quark mass m ^ 0. We regulate this 
by subtracting the trivial condensate contribution giving 

3 /■ , 2 • ' ""^ 



{(l(l)reg = J k (^siu 0fc - — j dk . (2.12) 

Our model is color confining and does not permit free solitary particles since the self- 
energy or dispersion relation 

2 f dq - ^ 

ek = msin0fc + kcoscpk - - J 'J2nf^(\^ " ^)(sin0fcsin0g + k ■ q cos 4>k cos 4>g) (2.13) 

obtained from the 1-body part of H, Eq. (p.l|) , is divergent (now there is no cancellation 
at the singular point k = q). Further, this divergence is also cancelled in the bound state 
equation but only for color singlet states (see below). Even though the self-energy diverges 
it is still useful to introduce the concept of an effective quasiparticle (constituent) mass, Mg, 
which can be extracted from the low momentum behavior of the gap angle. We introduce a 
running, dynamical mass, m(A;), by an effective Dirac spinor in canonical form 

Ulff{k)=N (2.14) 

E+M(k) -^^ 
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with normalization and E = \jM^{k) + k"^. Then using this equation and Eq. ( p.6| ) we 
equate the two relative normalizations between upper and lower spinor components yielding 
a relation between the running dynamical mass and gap angle 



l + sm(P{k) _ E + M{k) 

-, [Z.ID) 



^1 — sin (f){k) ^ 
or 

We identify the dressed quark or quasiparticle mass as = maxima M{k) and extract it 
from the low momentum behavior of the gap angle (essentially inverse of the slope near zero 
momentum). The value of characterizes the degree of chiral symmetry breaking and can 
be loosely regarded as the constituent quark mass associated with phenomenological quark 
models. 

Note that our expression for the running mass is functionally identical to the perturbative 
expression tana = ^. Related, since the rotated quasiparticle spinors have a running 
momentum dependence, they no longer rigorously provide a representation of the Lorentz 
group. Our form of dynamical chiral symmetry breaking violates Lorentz invariance which 
implies a preferred reference frame, namely the condensate rest frame. For most static 
observables such as masses, condensates and decay constants, Lorentz symmetry is not 
important. However, for some observables, such as electromagnetic form factors, care is 
necessary and boost corrections may be important. This issue is under investigation and 
will be reported in a future communication. 

C. Comparison to the Nambu-Jona-Lasinio model 



The classical effective model of Nambu-Jona-Lasinio (NJL) (see Ref. [jT3[ for review) 
entails various Lagrangian formulations, a common one being 

£ = ii|r + G[{^<i/){^^) - ('l'75^)(^75*)] (2.17) 



where G is a constant. It is customary to introduce the approximations ^'75\E' ^ and 
ilrxj/ ^ {^'^)vacuum, to huearize the equations of motion and then extract a constituent quark 
mass from the NJL mass gap equation. In this fashion chiral symmetry breaking is achieved. 
Our formulation extends beyond the NJL model in several important ways. 

1. Our approach is more general and permits explicit gluonic degrees of freedom (see Refs. 

The unification of the quark and glue sectors is crucial for a comprehensive 
treatment of hadron structure, especially for glueball and hybrid meson systems. 

2. Our formulation includes confinement and is renormalizable while the NJL model has 
neither. The NJL pointlike interaction would be recovered in the limit V{x — y) — > 
G6{x — y) which removes all important nonlocalities. 

3. Our model has a density- density interaction kernel with a different Lorentz structure, 
7o7o; which is the product of four-vector time components. As discussed in Ref. 

, a density-density (vector- vector) interaction is superior to the scalar-pseudoscalar 
displayed by the NJL model. 

4. The chiral symmetry breaking mode of the NJL is extremely restrictive yielding a 



constant quasiparticle mass, mdyn, and simple dispersion E = ^rrt^yn + k"^- Related, 
the NJL limit of our model also yields a more restricted gap angle since sin = 



Our method has a running mass and different quasiparticle dispersion 



which yields more realistic TDA and RPA hadron masses. 



III. NUMERICAL SOLUTION OF THE GAP EQUATION 

The gap equation ( p.9| ) has been previously solved for the harmonic oscillator potential, 
where it takes its simplest form as a differential equation (see Ref. [§,0), and also for the 
linear potential (see Ref. Here we summarize our analysis which confirms and extends 

the latter results. 
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To numerically treat the integrable IR singularity a regularization must still be imple- 
mented even though the final results are independent of this procedure. We considered 
several different regularizations. We first tried an analytical regulator (equivalent to a de- 
confining correction to the potential). This was an unstable algorithm and and conver- 
gence could not be achieved. We next examined the method of Ref. which off-sets the 
q-discretization by half a step in the kernel with respect to the k-discretization. This proce- 
dure was also rejected as it was less amenable for documenting the regulator sensitivity. We 
finally adopted the simplest method of omitting the singular point k = q. This also facili- 
tated a controlled senstivity study by just increasing the number of mesh points. Related, 
we adopted a variable mesh size to integrate more efficiently and mapped the integration 
variable g to f 

2 

^ Qmax 

q = 

for N points uniformly distributed in the interval v G (0, 1). 

Following Ref. we elected to solve the gap equation in form specified by Eq. ( p.lO| ) and 
also utilized the Gauss algorithm as described there. The Gauss method assures convergence 
but is rather inefficient for extensive sensivity studies in parameter space. We therefore 
modified our numerical approach by first finding a good approximate solution, ip^, to the 
non-linear gap equation and then obtained a linear equation for the desired correction, 6k, 
giving the final solution 



1pk = i'l + 



(3.1) 



to arbitrary accuracy. Substituting Eq. ( ^.1[ ) in the gap equation, Eq. ( p^.lOD , and dropping 
higher powers of 6k yields the approximate linear equation 



{k + 771)7/^^ - ^ /( 



oo q^dq 



J 2_ TOO q^dq Wr ( . fl ,/,0' ^ -^ ZLU '^fc^'. 



6. 



poo q' 



This equation is of the form 
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J dqA{k,q)5{q) = B{k) (3.2) 

which can be solved for 6k by matrix inversion. We found the Gauss algorithm sufficient for 
obtaining the initial approximate solution ip^. To achieve full convergence required up to 
12,000 mesh points, a factor of 60 more than the early calculations of Ref. 

We checked our computer codes by calculating two different toy kernels Vq = ^2(1^^2) ; 
Vi = and Vq = 0, Vi = ^2(1^^2) ; each designed to yield a known constant value for tp^- We 
then performed a series of cut-off sensitivity runs and mapped out the convergence rate as 
a function of mesh point number N which ranged from 100 to 12,000. We used the quark 
condensate as a test observable and also performed calculations for zero and non-zero current 
quark mass, m, with and without the Coulomb potential using as = .4. For m = 0, = 0, 
we determined the sensitivity to (the effective cut-off parameter) was slightly higher than 



previously reported @,|ri|] and given by 



1400\ 



113 —]MeV 



IT) 



3 



(3.3) 



Note that this number is somewhat smaller than the commonly accepted lattice value of 
about —(250 MeV)^. Including the Coulomb potential only increases the condensate to 
— (119 MeV)^. We therefore conclude an improved model ground state is needed which 
can be provided by including additional terms in the Hamiltonian, such as the quark-gluon 
minimal coupling (hyperfine) interaction. This point is also affirmed below in our RPA 
treatment which does yield a more realistic condensate value. 

Our other key result, which will be of interest in connection with chiral perturbation 



theory |14,TH], is for the constituent quark mass and the BCS condensate as a function of 



the u, d quark mass. Now it is necessary to use Eq. ( 2.12| ) and also impose an additional 



integration cut-off limit (g^ax around 10 GeV). This yields 

= 2.03 m 4-113.1 



Mg = 1.6 m-F 77.9 
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where the units are MeV. A precision calculation for Mg with m = yields the slightly 
higher value = 80.5 MeV (all values for the linear potential only). 

Finally, we note that our Hamiltonian is SUf{3) flavor symmetric, broken only by the 
small current flavored quark mass term. However, and quite significant, the vacuum proper- 
ties and gap angle exhibit substantial SUf{3) violations as evidenced by our strange quark 
calculations using a current mass of 150 MeV . Important violations occur even for a strange 
quark mass as low as 50 MeV. While this result may be model dependent it does suggest 
that certain chiral symmetry arguments in the literature regarding the strange quark sector 
should be taken with care. 

IV. MANY-BODY TECHNIQUES 

We now formulate mesons as excited states consisting of quasiparticles and seek approx- 
imate eigensolutions of our Hamiltonian. We first develop the TDA and then treat the RPA 
in the next subsection. Of the two, only the RPA preserves chiral symmetry, as we detail 
below. It is therefore more closely related to the Bethe-Salpeter formalism [|1^] incorporating 
the Schwinger-Dyson quark propagator using an instantaneous interaction in the rainbow 
approximation (equivalent to our gap equation). We will document this connection more 
formally in a future publication. 

A. TDA equation of motion 

The principle advantage of the TDA is that it is a controllable approximation which 
truncates the Fock-space expansion for a chosen level of calculational effort and resources. 
In terms of the quasiparticle operators introduced in Sec. II, we introduce the TDA meson 
creation operator 

QijpiTDA) = Y.j 7^^;r(^)5;(^ )Di{-k ) . (4.1) 
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A meson with quantum numbers nJP (radial-node number, n, total angular momentum, J, 
and parity, P) is then represented by the Fock space expansion 

I^™^)=QLp(™)I^) (4.2) 

containing a quasiparticle and quasihole excited from the BCS vacuum. The Hamiltonian 
equation is then projected onto this Ip-lh truncated Fock sector giving the TDA equation 

m'DA\[H.BiDl]\n) = {E^jp-Eo)Kf . (4.3) 

In evaluating the commutator we note 

[Ho,BiDl] = Q 

{n'DA\[H2,BiDl]\n) = {e, + e,)nf 
and for the two body potential 

■yS 

where is the BCS gap energy given by ( |2.13| ) and Hjy is the Hamiltonian component 
containing N field operators (after normal ordering with respect to the BCS vacuum). 

We can exploit the rotational invariance of our Hamiltonian and reduce the linear TDA 
equation to a one dimensional, nonlocal equation by an angular momentum decomposition. 
Introducing the orbital and spin angular momenta L and S, respectively, the meson state 
vector can be expanded in partial-waves involving a one-dimensional (radial) wavefunction 

*^r(^'')= E (^m^^m5|Jm,)(-l)^+^(i5i-7|^m5)irM^)^2f(fc) (4-4) 

LSniLms 

where again the color index is omitted. Note the phase factor and negative magnetic substate 
sign in the Clebsch-Gordan coefficient due to the transformation properties of antiparticles 
under the SU (2) rotation group. A thorough discussion is given in Ref. . 



Inserting Eq. ( [4.4|) into Eq. ([4.3|) yields the TDA partial-wave equation of motion 
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{Enjp - Eq- 2ek)'^l^s'{k){2J + 1) = J2AEmAmj,mjmLms{J'^j\L'^LSrns){J^rnA^rnT,\Jrnj) 

J dn,dn,Yr' mrm r ^^ni'^mnk - q\) 

where the function h'^g{k,q) contains the gap angle from contractions involving rotated 
spinors 

hjsif^^ l) = l[CkCg{SaSg'yl3 + S^l3gas) + (1 + Sfc)(l + Sq)das5p'y + (1 - Sg){l - Sk)ga69-tl3\ 

with 

Denoting the meson mass for state nJP by M^jp — E^jp — Eq and using the multipole ex- 
pansion formulas for the interaction yields the final TDA equation appropriate for numerical 
calculation 

(M„,p - 2ek)nf{k) = E r Kiis^ik, q)nif{q)^, ■ (4-5) 

AS -^0 127r^ 

Note the Hamiltonian spin dependence generates a kernel that couples different orbital and 
spin states. 

We now apply these equations to the low lying meson spectrum with quantum states 
specified by I'^{J^'^) having C parity, C = (-1)^+'^, and G parity, G = (-1)^+'^+''^. In our 
model we neglect the small electromagnetic (isospin violating) effects as well as coupling to 
the gluon sector so that / = and 1 states are degenerate for the same J^^ . For pseudoscalar 
states, J^*^ = O""*", S — L — J — giving only one wavefunction component (no coupling). 
This is also the case for scalar mesons (L = 1,5" = 1) having J-^^ = However for 

the vector meson sector J^*^ = 1 both {L = 0, S = 1) and {L = 2, S = 1) waves are 
allowed and, in general, will be coupled. Similarly for low lying pseudovector mesons having 
(L = 1, 5" = or 5" = 1) and tensor mesons, with (L — 1,S — 2) and {L — 3,S — 2), 
there will be coupled equations. Although these equations are not difficult to solve, we 
only include the lowest orbital partial-wave component and neglect all coupling since it has 
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been computed small for the harmonic oscillator potential . There is then only one kernel 
K{k, q) for each meson state with structure: 

• pseudoscalar, L = S* = J = 0, 

K{k, q) = 2{ckCgVi + (1 + SkSg)Vo) 

• scalar, L = S = 1, J = 0, 

K{k, q) = 2{ckCqVo + (1 + SkSq)Vi) 

• vector, L = 0, S=J=1 (neglecting the tensor L = 2,S = J = 1 coupling), 

K{k, q) = 2c^CqV^ + (1 + su){l + s,)^ + (1 - - Sk){ ° ) 

• pseudovector, L = J=1, 5' = (degenerate with L = S = J = 1), 

K{k, q) = cucq{% + V2) + 2(1 + SkSq)V^ 

• tensor, L = S = \, J = 2 (neglecting L = 3,S=J = 2 coupling), 

K{k, q) = CuC,{W2 - Vo) + (1 + Sk){l + Sq)Vr + (1 - Sk){l - ,j ^^^3-7Vi 

Note for the pseudovector mesons the kernels for 5* = and S = 1 are identical which differs 
from Ref. §. 

We have also applied our approach to other flavored (s and c) meson systems and have 
obtained similar, but more complicated TDA equations. As a representative result, consider 
the pseudoscalar D meson with a u (or d) and c quark. The gap equation for the c quark 
remains the same, except for current mass, now 1.2 GeV ^ which gives a different gap energy, 
e^, and angle, s^. The TDA equation, however, has a different form and generalizes to 

{Mn -el- el)^o{k) = \ j ^.^niq)- 



VT+4xn + s^VT+4xn + s'. + J I - s"v^r^v/r^,/l - V^{k, q) 



1 - s^v/r^v/TT^Vl + sl + v^TT^y'l + s^v^r^y'l - si ) V^{k, q) 
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with obvious form for other mixed flavors. All equations are finite for k = q as the IR 
divergence terms from the confining potential again cancel. We have also derived and solved 
the TDA equations for other spin parity states which is further detailed in Sec. V. 

B. RPA and the quasiboson approximation 

The TDA can be improved by utilizing a better vacuum with additional quasiparticle cor- 
relations beyond the BCS. Consistent with many-body applications in other disciplines we 
now formulate the RPA |]TB|,|TP[ and introduce a new vacuum, \^Irpa), having both fermion 
(two quasiparticles or Cooper pairs) and boson (four quasiparticles or meson pairs) correla- 
tions. The RPA meson state 

mi'l) = Qljp{RPA)\QnPA) (4.6) 
involves a meson creation operator which is a generalization of Eq. ( [4. 1|) 




The RPA vacuum then satisfies 

QnjpiRPA)\QRPA) = 

which, because of additional correlations from admixtures of particle-hole excitation states, 
is not true for the BCS vacuum. 

To derive the RPA equations of motion we use Eq. ( |4.3| ) and replace the BCS vac- 
uum with IQrpa) and also substitute ^E'^p^ for ^E'^zm generate one equation for the X 
component. We then repeat, changing the commutator to [H, BaDp] to obtain the Y equa- 
tion. Following standard treatments in other fields of physics, we also invoke the quasiboson 
approximation and treat the fermion pair operator BD as a pure boson operator. This sig- 
nificantly reduces the commutator algebra complexity and generates one of the two coupled 
equations for the RPA wavefunctions X and Y . 

For the important pseudoscalar meson channel we obtain for the excited state n 
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2e,X^{k) + 1 ^ [X'\q)F{k, q) + Y-{q)G{k, q)] = M^X^\k) (4.8) 

2e,y"(fc) + 1 ^ [Y^iq)F{k, q) + X^{q)G{k, q)\ = -M^Y^{k) (4.9) 

where 

F{k, q) = 2c5CfeVl + 2(1 + SqSk% (4.10) 

G{k, q) = 2cgCkVi - 2(1 - SgSk)Vo ■ (4.11) 

Similar expressions directly follow for the other spin-parity states. 
We adopt the standard normalization for the RPA wavefunctions 

= {nRPA\Qu'Ql\nRPA) = 5,u' (4.12) 

yielding 

/ ^dkiX" ikYX^ik) - ikyVik)) = (27r)2(5,,, . 

^0 

The RPA equations, which reduce to the TDA equations in the limit y or G ^ 0, are 
again an eigenvalue problem for M^jp which can be easily diagonalized. Related, the matrix 
size can be reduced by a factor of 2 using the variables X -\-Y and X — Y. Finally, the 
equations are also IR finite for the sigular point k — q. 

C. Weak decay constants 

A crucial test of any approach is the ability to describe hadronic decays. In this paper we 
compute weak decays and defer our analysis of hadronic decays to a subsequent publication. 
For a pseudoscalar meson P with momentum p^, energy Ep, mass Mp, the weak decay 
constant, /p, is defined for our normalization by 

(O |A^(0)|P(p )) = -^fpp, . (4.13) 

Here A^^x) = ^(x)7^75^(x) is the axial current which specifies the chiral charge operator 

Qb = j dxAo^x) = j df^^(f)75^(f) . (4.14) 



Simplifying Eq. ( [4.13| ) for a meson at rest yields 



fp = -^m^oh.^mpm . (4.15) 

Applying this result for the TDA pion wavefunction gives the TDA pion decay constant 

1 



/™^ = ^=^— / nDAiQ)s, q'dq . (4.16) 

Similarly, for the RPA pion wavefunction we obtain 

1 

= / s,{X-{q) + Y-{q)) q'dq . (4.17) 



7r^(27r)3M^ 

Our results easily generalize to the SUf{3) flavor nonet. Now there are nine axial charges 
given by 

Qt = J dxA^,{x) = J »t(f)^^^vi>(f) 



where the eight Gell-Mann A" matrices are supplemented by A° = y |/ to obtain both the 
octet and the singlet under SUf{3) transformations. The appropriate generalizations of Eq. 
( [4.16| ) are then: 



27r^(27r)3M/^ 

1 r'^ 
37r^(27r)3M^3 

37r^(27r)3M,„ 

We will use the above results in the next subsection to derive a generalized Gell-Mann- 
Oakes-Renner relation and also in Sec. V were we report numerical results. 

D. Chiral symmetry and the Gell-Mann-Oakes-Renner relation 

Chiral symmetry is a significant element of hadronic QCD and should be present in 
all realistic models. Even though our vacuum properly exhibits dynamic chiral symmetry 
breaking, our model Hamiltonian does indeed respect this symmetry since the commutator 



[H,Q^] = J (if\&"^(f)\&(f) ^ 
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essentially vanishes, consistent with the small u, d quark mass. Related, our RPA states 
also preserve chiral symmetry as the RPA meson creation operator commutes with the chiral 
charge in the chiral limit (m 0) 

[Q\RPA),Q,] = Q . 

However, the TDA operator, Eq. ([4.1D above, does not commute with and violates chiral 



symmetry since it is not fully symmetric in operator structure (only contains B'^ D'^ and not 
BD). This can also be documented by chiral transforming the TDA meson state verifing 
that B'^ D'^ rotates to combinations of B'^B, DD'^ and DB. Hence, the TDA ansatz is not 
closed under a chiral rotation and Goldstone bosons will not appear in the TDA spectrum. 
We therefore expect significant, but unphysical, chiral symmetry violations in the TDA 
calculations and anticipate the TDA pion mass to be much larger than in the RPA which is 
confirmed in Sec. V as only the RPA calculations yield a Goldstone pion in the chiral limit. 

On the other hand, the RPA ansatz, Eq. ( [4.7| ), is chirally invariant since it is closed 
under this rotation in the chiral limit = Hence the mechanism of spontaneous 

symmetry breaking is present and the pion mass is zero according to its Goldstone boson 
nature as we numerically verify in the next section. 

With these results we now derive two different Gell-Mann-Oakes-Renner (GMOR) re- 
lations, one based upon exact model eigenstates while the other relates RPA states and 
the BCS vacuum. Assume we have the complete set of exact eigenstates, \n), to our QCD 
model Hamiltonian, including the vacuum ground state Inexact)- Evaluating the double 
commutator 

{^exact\ [Qb, [Qb, H]] \Vtexact) = 4:m{qq) exact (4-18) 

then generates the exact quark condensate. Evaluating the double commutator again, but 
now invoking twice the completeness relation, 1 = J2n\''^){My ^^^1 identifying the decay 
constant relation, Eq. ( [4.15| ), leads to the generalized GMOR relation 

-2m{qq) = Y.Kfl (4.19) 

n 
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summed over all (ground and excited) pseudoscalar meson states with mass M„ and decay 
constants /„. 

This can be extended to flavor using two of the nine SUf{3) axial charge operators, 
and Qg, to derive the following GMOR relations: 

• a= b = 1, 2, or 3 

TT" 

• a = b = 4, 5, 6, or 7 

\ ^ J j^n 

• a = b = 8, 



• a = 




Because the exact eigenstates will generally not be available, these equations are of lim- 
ited value. However, they do provide testing criterion for approximation solutions. Further, 
since decay constants are suppressed for excited states n we can drop higher terms to obtain 
more useful relations involving ratios such as 

MKfiA ^ _ rriu + nis ( {uu + ss) \ 
MJ^ J ~ 2m, 1^ 2{uu) J 

which we will use in Sec. V in connection with our discussion of the kaon mass calculation. 
Next we utilize Thouless' theorem applied to the chiral charge operator 

{Q\ [Qs, [Q5,H]] \Q) = 2Y,\{^^^'p^^\Q,\nnpA)\\E^ - Eo)rpa (4.20) 

n 

to immediately derive the RPA GMOR relation 

- 2m{qq)Bcs = {M'J^)rpa (4.21) 
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where we have again dropped the excited state decay constants. Note that the left hand side 
entails the BCS vacuum while the right hand side involves the RPA states and energies. This 
relation clearly predicts the RPA pion is a Golstone boson in the chiral limit (i.e M^r 
for m ^ 0) which we numerically confirm in the next section. 

Finally, calculating the RPA meson mass spectrum does not require obtaining the RPA 
vacuum (see Eq. (|4^ )), but computing the RPA decay constants does. Determining IQrpa) 
is actually quite difficult, however, the leading correction can be approximately calculated 
using another theorem by Thouless [T^ 



I^rpa) ^ \n) + {F^y\^) (4.22) 
where is a TDA type operator given by 

a/3 '' 

with the fermion pair operator B^D^ now obeying bosonic commutation relations. Here /* is 
an unknown amplitude assumed to be small (this is not true in the chiral limit as discussed in 
the next section). The approximate RPA vacuum is thus described as a mixture of the BCS 
vacuum having Cooper pairs (Eq. p.7| ) and two quasibosons (mesons) coupled to vacuum 
quantum numbers (0^^). Since, as shown in the next section, the RPA is of importance 
for mainly the low-lying pseudoscalar states, we assume that only these states contribute 
to \^rpa) and neglect all others. Our result easily generalizes to scalar or other mesons 
(appropriately coupled to = 0++). Within this approximation we obtain an improved 
quark condensate to be compared with Eq. (|2.11| ) 

mlRPA ^ m)BCS + 



i + 2Ea/3/rfg1/a/3(g)|2 ■ 

To obtain /q^, we impose Qt,{RPA)\VLrpa} = and use Eq. ([4.22|) neglecting higher 



order terms corresponding to Fock states with more than two pions. We also only retain 
ground state meson (pion) contributions from X and Y yielding 
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Uiq) = iV-iF(g)(-l)^+"+'3(l - a^m) (4-23) 



with a normalization constant depending on both RPA wavefunction components 



The improved RPA quark condensate is 

{qq)RPA ^ {qq)BCS + c / SqY{qfq^dq {AM) 

Jo 

where the constant c is given by 

2JY{q)Ydq 



UY{q)X{q)qHqf + \Uy{qYq^dqY ' 
Also the improved pion decay constant to this vacuum is 

U ^ I ^ / s,{X{q) + c7i^N^Y{q))qHq (4.25) 

7r^(27r)3M^ J 

in contrast to Eq. ( [4.17|) . We have found that for meson masses above 800 MeV Y is very 
small and there is no essential difference between RPA and TDA. 

We now comprehensively apply the above formulas and conduct a comparative analysis 
of the TDA and RPA approaches. 



V. APPLICATIONS AND NUMERICAL RESULTS 

In this section we first present and discuss our TDA spectra for the pseudoscalar and 
vector mesons and then in subsection B we compare to our RPA results for both meson 
masses and the decay constants. Subsection C treats mixing for the t] and rj' mesons while 
subsection D details applications to charmed mesons, especially J/\E' states. 



A. Tamm-DancofF spectrum 



Solving the TDA eigenvalue problem is straight forward. The diagonal part of the matrix 
contains the IR singularity that is rigorously controlled by cancellation, again permitting 
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numerical regulation by simply skipping the point g = as in the gap equation solution. 
Due to the linear nature of this system, results are convergent for a mesh as sparce as 700 
points. 

Tables I and II summarize the resulting TDA meson spectrum corresponding to the five 
JLS kernels specified in Sec. IV A. The energy difference between pseudoscalar and vector 
states is about 200 MeV for the u/d quark mesons, 80 MeV for the open flavored {K, 
K*) and only 50 MeV for the pure strange composites (0). Our one parameter model can 
not accurately describe the entire observed splitting indicating additional dynamics beyond 
simple spin interactions from a Dirac spinor field is needed. In general, the masses are in 
good agreement with the PDG accepted values for the vector mesons, but the TDA low- 
lying pseudoscalar sector is deficient. This is expected since in this channel vacuum (chiral 
symmetry) effects are most prominent. In the scalar channel the situation is more confusing, 
since other hadron states, some with explicit gluonic structure, can more easily mix. Since 
our unified model allows us to treat glueballs, mesons and hybrids comprehensively, future 
work will further address understanding this channel. Interestingly, our lightest /o mass, 
which has a P wave oribital excitation, is below 1 GeV . We also find that the computed 
TDA masses are not very sensitive to details of the vacuum that enters via the gap function 
characterizing the BCS ground state. 

As described above, we may also use a simplified TDA equation to extract a constituent 
quark mass. In the chiral limit, the mass obtained is 51 MeV and this dressing is roughly 
constant, consistent with SUf{3) symmetry, up to current masses of 150 MeV, where the 
generated constituent mass is 203 MeV. We are therefore dealing with light dressed quarks, 
even for strange quarks which may prove a deficiency when calculating electromagnetic form 
factors. We shall also address this issue in the near future. 

Another important feature of our relativistic effective Hamiltonian is that it naturally 
includes the kinematical and spin dependent interactions (e.g. spin-spin, spin-orbit, tensor). 
These effects are very important in the light quark sector even after chiral symmetry break- 
ing, because the light quark constituent mass generated in our scheme, around 80 MeV, is 
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still small as compared to our interaction scale (424 MeV). In particular, notice in Table II 
the spin splitting between the 1^+ and 2++ mesons, all having the same L and S quan- 
tum numbers (also observe the large radial excitation in each channel). The level spacing is 
consistent with Ref. but very different from the naive expectations from the constituent 
quark model. The is due to the confining (nonperturbative) L ■ S coupling, which is only 
of order as in the constituent quark model. The spin spacing is governed by the matrix 
element 

{LSJ\L ■ S\LSJ) = ^[J{J + 1) - L{L + 1) - SiS + 1)] 

which, for L = S = 1, reduces to —2 + ^^^^^^1) . The —2 in this expression describes the light 
scalar meson, whereas the J{J + 1) contribution explains the splittings in Table II. The 
calculated 2++ mass is much heavier than the lightest observed /2 at 1270 MeV. Clearly, 
including coupling to the L = 3, S = 1 channel as well as multi-quark Fock states will alter, 
but improve, our prediction. Further, and as mentioned previously, additional quark-gluon 
interaction terms should also be included in our Hamiltonian which have a different Lorentz 
spin structure. In our model this would generate a weaker hyperfine interaction of order as- 
This would also provide additional splitting between 5* = and S = 1 levels that would 
improve our TDA, and especially RPA, vr-p mass underprediction. Such an analysis would 
fully clarify the relative importance of chiral symmetry versus spin interactions as the latter 
is generally attributed the dominant effect in conventional constituent quark models having 
color magnetic, effective one gluon exchange potentials (see Ref. [0). We are currently 
examining this issue and will report results in a future paper. 

Finally and also related is the spin-orbit splitting for other flavors which is summarized 
in Tables I and IV. In Table IV we illustrate the familiar J{J + 1) dependence by fitting the 
TDA spectrum for different flavors. 
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B. RPA spectrum and decay constants 



We now present our RPA results and compare with both TDA and observation. Table 
III and Figs. 2 through 7 highlight our key results. In general the RPA and TDA masses 
agree except for the light pseudoscalar mesons where the RPA provides a better description. 
This is because, as discussed above, only the RPA correctly implements chiral symmetry 
as illustrated in Fig. 3 where the lightest scalar and pseudoscalar masses are plotted as 
a continuous function of the current quark mass. Note that only the pseudoscalar mass 
approaches zero in the chiral limit consistent with Goldstone's theorem (the observed pion 
mass is reproduced for a u quark mass of about 2 MeV). 

The kaon system reveals the largest model deficiency (see Table I). Even the RPA in 
the chiral limit produces a too massive kaon, about 850 MeV. As indicated by Fig. 4, to 
reproduce the observed kaon requires a strange quark mass of about 50 MeV. A detailed 
analysis reveals that the explicit contribution of the current quark mass to the RPA equation 
-through its appearance in e{k)- is only additive. It is the gap angle, introducing an implicit 
flavor dependence, which inhibits a lighter kaon mass. We could fit both pion and kaon 
masses by adjusting the current quark masses, but we prefer awaiting improvements from 
renormalizing the quark gap equation (work in progress and will be subsequently publish). 
The gap angle for a non-zero current quark mass is very sensitive at high momentum to the 
dominanting Coulomb potential and sizeable corrections are expected. Also for higher lying 
excited meson states there is also the issue of two particle, two hole Fock state contributions. 

In Fig. 5 we compare our scalar, pseudovector and tensor meson TDA and RPA predic- 
tions to data. In general there is qualitative agreement. Note that for the higher excited 
states above 1 GeV the TDA and RPA results are identical and it is clear that these systems 
are not governed by chiral symmetry. 

We also illustrate the behavior of various wavef unctions. In Fig. 6 we detail the difference 
between TDA wavef unctions for the /o and tt mesons. Figure 7 compares the TDA and RPA 
wavef unctions for the pion. 
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Finally we discuss decay constants. Both the TDA and RPA pion decay constants are too 
small, about 17 MeV, in contrast to the observed value of 93 MeV. This is consistent with 
the Orsay results. We again attribute this to the model Hamiltonian and vacuum as the 
BCS angle does not have sufficiently high momentum components. However, if we use the 
approximate RPA vacuum in the quasiboson approximation (see Eq. (4.25)) the improved 
RPA decay constant increases to 57 MeV. Appropriately, the condensate also significantly 
increases to —(320 Me \^)^, in much better agreement with lattice results ~ —{250MeV)^. 
The latter result is consistent with applications in nuclear physics where the RPA tends to 
overcorrelate the ground state. 

Since we have only approximately evaluated the RPA vacuum by truncating at the two 
pion level it is not surprising that does not agree with measurement (or chiral perturbation 
theory results) and therefore needs further refining. We expect the truncation to the two pion 
level, Eq. ( |4.22|) , to be reasonable provided we are not in the chiral limit, since then Y <^ X. 
However, we would like to be able to calculate in that regime. Further, the calculation in Ref. 
PB| points to a necessary decrease of the BCS condensate when including coupled channels. 
They argue that for a pure chiral pion, coupling additional channels would decrease its mass 
which might even become negative, destabilizing the vacuum. This argument seems sound 
and since we are above the chiral quark mass limit, at a model value m = 5 MeV where our 
pion is too massive (277 MeV), this decrease is a welcome improvement. We defer further 
discussion until publication by our collaborative effort with the Lisbon group which will also 
clarify pionic correlations in the ground state of our approach. 

Recalling that our computed condensates require renormalization except in the chiral 
limit, we can only test the GMOR relation for m = which is trivially satisfied in the 
RPA. However, we note that the decay constants for excited pion states are much smaller 
and rapidly approach zero in the chiral limit. Hence the GMOR relation is satisfied with 
predominatly the first state. This is not true for heavier quarks and to numerically satisfy 
the GMOR requires including several eigenstates. 
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C. r] — rj' mixing 



An extremely challenging but still not understood problem is the t], rj system and at- 
tending flavor mixing of the light quarks. Although our effective Hamiltonian has an explicit 
flavor dependence through the current quark masses, it still conserves flavor. However, if 
the gluon sector is included, such as through the hyperfine, minimal coupling interaction, an 
effective flavor dependence naturally emerges though higher order quark-glue-quark effects 
and dynamic mixing of flavor states is possible. We are currently deriving such a term which 
is similar, but more rigorous than the t'Hooft interaction based upon instantons (classical 
glue). This will be reported in a future communication, however, it is still of interest to 
perform a simple 77, rj mixing analysis by introducing a flavor off-diagonal interaction as we 
now detail. 

With no dynamic flavor mixing, the t] and Vj are (poorly) modeled as SUf{3) octet, r^g. 



and singlet, rjQ, states respectively given by (we adopt the convention used in Refs. p^ , |25| 



V^Vs = ]j^nn- ^j'^ ss = cosOsUfis) nn - sinOsUfis) ss (5.1) 

v' ^ Vo = nn + ^ ss = sinOsUfis) nfi + cos9sUf(3) ss (5.2) 

involving the isoscalar nn = [uu + dd)/\/2 and ss states and pure SUf{3) mixing angle 
^sUf{3) = 54.74°. We have computed the TDA masses of the pure nn and ss meson states to 
be Mn = 612 MeV and M, = 1002 MeV, respectively. Hence the predicted, pure SUf{3), 
1], 1]' masses are 

^ = cos^ esuf{3) M„ + sin^ OsUfis) = 872MeV (5.3) 
M^, ^ = sin^ ^5^^.(3) M„ + cos^ Osu^is) = 7A2MeV. (5.4) 

For the RPA, M„ = 290 MeV and M, = 978 MeV, yielding M^ = 749 MeV and M^. = 519 
MeV. Similarly the SUf{3) t], rj' decay constants are given in terms of the u/d isoscalar, 
/„, and strange quark, /g, decay constants 
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fr, ~ fm = cos^ 0suf{3) fn + sm^ 9sUf{3) fs (5.5) 
/r?' ~ fvo = Sin^ OsUf{3) fn + COS^ OsUf{3) fs- (5.6) 

Note the quadratic dependence on angles due to expressing both the axial current and meson 
states in their SUf{3) representations. Using the TDA computed values of = 17 MeV 
and fs = 75 MeV, yields /^g = 56 MeV and /^^ = 40 MeV. The RPA values are /„ = 57 
MeV and = 75 MeV, giving ^ = 68 MeV and ^ = 63 MeV. 

We can now improve these results by generalizing our Hamiltonian, which is diagonal 
in flavor space, to have off-diagonal matrix elements. The simplest prescription is to just 
add a constant, {qq\H'\qq) = A, to both diagonal, {nn\H'\nn) = 2A, {ss\H'\ss) = A, and 
off-diagonal {nn\H'\ss) = \/2A, terms giving 

^M„ + 2A V2X ^ 
Diagonalizing H leads to the new, mixed mass eigenvalues 



(5.7) 



M„ + M, + 3A _ 1 + + 9A2 + 2(AM„ - M,M„ - M,A) (5.8) 
M,, = ^n + M, + 3A ^ 1 + ^2 + 9p + 2(AM„ - M,M„ - M,A) (5.9) 

and eigenstates 

r] = cos{9sUf{3) + Op) nn - sm{9suf(3) + Op) ss = cos9p r]s - sin 9p rjo (5.10) 
T]' = sin(6l5[/^(3) 6lp) nn + cos{9sUf{3) + 9p) ss = cos 9p rjo + sin 9p rjs (5.11) 

involving rotation by an additional angle 9p that is a function of A. The mixed, presumably 
more physical, decay constants are then 

= cos 9p /^g - sin 9p fr,o (5.12) 
/^/ = cos 9p + sin 9p /^g . (5.13) 

Performing a least squares fit to the observed masses (M^ = 547 MeV , M,f = 958 MeV), 
yields A = -33Me\/ (^p = -61°) for the TDA, which in turn produces = 541 Me\/, 
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= 974 MeV, fr, = 62 MeV and = -30 MeV. For the RPA, A = 82 MeV {9p = 
-44°) generating = 433 MeV, M^, = 1081 MeV, fr, = 93 MeV and f^, = -3 Me\/. 
It is interesting that while the simple mixing provides improvement, the TDA masses are 
in better agreement than the RPA. We could also improve the decay constants utilizing a 
two angle mixing formalism but refrain since clearly a more sophisticated treatment is 
necessary which will be provided by our quark-gluon coupling formulation in the near future. 



D. Heavy mesons 

The constituent quark models and non-relativistic expansions of QCD offer more reliable 
results for heavy quark systems where physical intuition from Quantum Mechanical bound 
states is more appropriate. Hence the charmed mesons afford a good limiting testing for 
our relativistic approach. We again calculate the spectrum of the charmed mesons (see 
Figures [§,0) and of charmonium (Fig. |Ty]) using our many-body model. The TDA is now 
sufficient since chiral symmetry, crucial for the light mesons, is not a constraint and the 
RPA will produce the same results. Using a charmed quark mass of 1200 MeV, the general 
features of the spectra are well reproduced and the radial excitations, \E'(2S') and rjc, are 
adequately described. We therefore expect our predictions for the remaining unconfirmed 
states to be reasonable. 

The angular momentum splittings of these systems are known to be dominated by the 
one gluon exchange potential (OGE) which we have not included in this calculation. Hence 
there will be improvement from future calculations based upon our renormalized project. 
A general feature reflected by our charmed spectra is the near vanishing of the spin-spin 
interaction, leading to degenerate 0""*" and 1 states. Since the 100-150 MeV hyperfine 
splittings in these systems will presumably be recovered when we include the perturbative 
OGE, we do not comment further. As for the spin-orbit splittings, our results are too large 
for the D mesons and too small for the Xc mesons. These splittings are adequately explained 
in non-relativistic quark models (see Ref. |^) where the absence of a large spin-orbit effect 
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for light quark masses is attributed to the cancellation between the Thomas precession in 
the confining potential and the one gluon exchange effective potential, although it is a bit 
concerning that the actual splitting between the Xco and the Xci is twice the size of the 
Xci-Xc2 splitting, when according to the spin-orbit J{J + 1) rule, it should be a half. 

Finally, we note that by including D waves in the charmonium spectrum we are able to 
resolve the long standing "overpopulation" problem of J/^ states relative to observation. 
This is clearly illustrated in Fig. 10. This (previous) deficiency in number of cc states has 
been characterized as evidence for glueballs and/or hybrid mesons because the J/^ system 
is believed to be gluon rich. Our result suggests, however, that simple level counting may 
not be effective in identifying hadrons with explicit gluonic degrees of freedom. 

VI. OUTLOOK AND SUMMARY 

Before summarizing our results we comment on the strengths and weaknesses of our 
many-body approach as well as some attending, open hadronic physics issues. Beginning 
with our Hamiltonian, H, the current-current (density-density) color interaction forbids free, 
isolated colored objects in the theory. For the vacuum this is realized in the BCS by an 

A — * 

infinite shift in the free quark self-energy due to the integral of V{\k — q\). Similarly for 
hadrons in both TDA and RPA, colored composite objects (e.g. diquarks) are precluded by 
the appearance of V'(O) which is divergent, whereas in the singlet channel this divergence 
is removed by vanishing color factors. Next we note that H conserves chiral symmetry yet 
our BCS vacuum properly exhibits dynamic chiral symmetry breaking. Further, our RPA 
pion emerges as a Goldstone boson in the chiral limit. This is not true for the TDA pion 
since only the RPA excitation operator commutes with the chiral charge. Another, signifi- 
cant model feature is that Fock state truncation is a controlable approximation amendable 
to systematic improvement. Thus our Hamiltonian many-body approach is an attractive, 
promising method for comprehensively investigating hadronic structure as it embodies con- 
finement, chiral symmetry breaking and orderly construction of multiparticle Fock states. 
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It also provides an excellent vehicle for testing more fundamental effective Hamiltonians as 
well as affording a powerful phenomenological framework for hadron structure. Considering 
the form of if, with only a single predetermined dynamical parameter, it is encouraging 
that the chiral limit is adequately reached in the RPA and that the meson spectrum is in 
qualitative agreement with experiment. To achieve detailed, quantitative descriptions will 
require further improvements in both the Hamiltonian and effects from including higher Fock 
space components. In particular, both the high energy behavior (Coulomb potential) and 
quark-gluon coupling effects (efectively instantons) will be incorporated and reported in a 
future publication. Finally, our current study is similar, but more extensive than the Orsay 
analysis due to our application to multifiavor systems. Our results are also more realistic 
(and numerically more difficult) then that work since we have utilized a linear confining 
interaction, determined by lattice and Regge phenomenology, which generates complicated 
nonlocal integral equations, rather than solving a simpler differential equation for a harmonic 
oscillator potential. 

Summarizing, we have performed approximate, but large-scale diagonalizations of an 
effective Coulomb gauge Hamiltonian utilizing standard many-body techniques. Using the 
BCS, a non-linear gap equation has been derived and accurately solved to provide vacuum 
properties (quasiparticles and condensates). Incorporating only predetermined parameters 
(string tension, a, and reasonable current quark masses), we have qualitatively reproduced 
the low energy u, d, s and c meson spectra. Most importantly, we have obtained a chiral 
pion, detailed that chiral symmetry is responsible for the large vr — p mass splitting and 
resolved the problem of overpopulation of theoretical J/\l/ states. 

Future work will address the full Hamiltonian in the combined quark and gluon sectors. 
In particular, we will obtain improved spin (hyperfine) and flavor (t'Hooft) interactions from 
quark-glue coupling. This should provide a better description of vacuum properties and the 
scalar /pseudoscalar masses, especially the vr, K, rj and rj'. We will also include more complex 
2 quasiparticle-2 quasihole Fock states for heavier u/d mesons as well as 3 quasiparticles for 
baryons and hybrids. Much of this work is in progress and will soon be reported. 
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APPENDIX A: BCS VACUUM STATE 



We further discuss the relation between the BCS rotated, and the trivial or perturba- 
tive, |0), vacua. We first note that the BCS vacuum state given by Eq. ( p.7|) is not a unitary- 
transformation and does not have a finite normalization. This is because the operator in the 
exponential is not antihermitian. It is therefore necessary to normalize matrix elements by 
dividing with and this is implicit in our presentation. Alternatively, and equivalently, 

\Q) can be represented by a norm preserving unitary transformation of the form 

\Q) = e^'-^l 0) 

where 

A^=T. J dktandkM^,^,bl{k)dl{-k) 
A1A2 

and all flavor and color indices are suppressed. Here Ma^a2 cire matrix elements of the Pauli 
matrices M = a ■ k 

It is interesting to note that the BCS vacuum state is orthogonal to the trivial vacuum, 
{0\Q) = 0, in the infinite volume limit. Further, the Hilbert space vectors constructed from 
the two different vacua are also orthogonal provided the BCS angle is nonzero (vacuum 
condensates are present). Because of this property the BCS rotation has been called a 
pseudounitary transformation 0. Consult this reference for further details (note they have 
a different phase convention and use M = a ■ k{ia2))- 
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APPENDIX B: GENERAL TDA EQUATION 



Using the phase convention of Ref. |28], the general TDA meson equation for arbitrary 
angular momentum is 



Jo 



q^dq 



vr^ 



;i + Sfc)(l + s, 



"vK{k,qWd!'{q)+ 



12 



ILSf 1 lLSfghLiL2 2 _ 

where the angular momentum products are 



n = (10L0|/0)(10A0|/0)(2/ + 1)V(2S + 1){2S + 1)(2A + 1)(2L + 1) 
1 

■W{E^S^; i/)W^(LlAl; //)H^(L^AS; J/)(-l)^+^+i(l + (-1)^+^) 



and 



n = (2/ + l)(2^ + l)(2/i + 1) V(2Li + 1)(2L2 + 1)(2L + 1)(2A + 1)(2S + 1){2S + 1) 

2 

.(_l)(i+2(/-^)+i-A+i-J) ^1010|LiO) (IOIOIL2O) (LiOAO|0) (L2OLOI/O) 

lS)iy(L2lii; 

•W^(L2LLiA; Z/i)l^(SLSA; Jh)W{fi:gS; h)W{L2gLJ; ^h). 

These formulas have been applied to several different meson spin, parity states as presented 
in Sec. IV. 

The moments of the angular integrations for the linear potential are obtained from 

where x = k ■ q. These can be calculated using the recurrence relation 

k'-r Ana d I I 



" q q dk [ 2qk li a; - ^ 

or by explicit integration 
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" q 2g2 dk 



k \ n — l n — 2 



+ ■ ■ ■ + XUJ 



n-2 



'CO -V 



k2 I 2 

where u = . Evaluating the first three moments yields 



yL — IGtto" 



y L 27r(T 



V2 = Sttct . ^ „ In 



k^q^ \k — q J k'^q^ {kp' — q^Y 
For the Coulomb potential, Vc — ]^^) ^ — ' angular integrals are 



-27rQ;5 , I ( k + q' 
m 



k — q 



''I — — ; — I — ~i — "^o 



qk 2qk 



2qk 
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TABLES 



TABLE I. TDA ground and first excited states in MeV. Linear potential only, a = 0.18 GeF^. 





0-+ 


0++ 


1— 


l+± 


mi 


m2 


El 


E2 


El 


E2 


El 


E2 


El 


E2 








586 


1473 


817 


1667 


798 


1602 


1076 


1818 


5 


5 


612 


1494 


850 


1675 


800 


1615 


1093 


1835 


5 


10 


624 


1503 


861 


1703 


803 


1619 


1100 


1843 


5 


150 


877 


1679 


1086 


1873 


957 


1743 


1273 


1988 


150 


150 


1002 


1808 


1297 


2044 


1044 


1849 


1416 


2116 



TABLE II. TDA u/d mesons (L = 5 = 1, no L = 3) in the chiral limit. 



jnC 


El 


E2 


E2, 


0++ 


817 


1667 


2301 


1++ 


1076 


1818 


2411 


2++ 


1767 


2281 


2749 



TABLE III. ( •liirMl H\'irimc(r}' V^rc^Mkiiig in tlic; RPA: scMlrir \-s. pscudoscalar tij")tx-(ruiri. 



mi 


m2 
















729 


1652 





1435 


5 


5 


775 


1679 


300 


1463 


5 


10 


794 


1641 


350 


1475 


14 


14 


838 


1719 


441 


1502 


150 


150 


1288 


2042 


978 


1805 
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TABLE IV. TDA fits to the spin-orbit splitting, Ej+i - Ej = AJ{J + 1), for L = S = 1 as a 
function of the current quark mass. 



mi/m2 


5/5 


150/150 


5/1200 


150/1200 


1200/1200 


A {MeV) 


162 


77 


74 


40 


10 



40 
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FIG. 1. Quasiparticle energies for the u/d,s quarks and the gluon [||. 
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FIG. 2. Pseudoscalar and vector(L = 0) TDA and RPA meson spectra. Only the RPA provides 
the correct chiral limit for the pion. 
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FIG. 3. Chiral symmetry in the RPA. For m — >^ the pseudoscalar (sohd) but not scalar 
(dotted) meson mass vanishes. 
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FIG. 4. Chiral behavior of the kaon and first radial excitation. Light quark mass is MeV. 



44 



E(MeVY 



2000 



1600 



1200 



800 



400 



f (2060) 
n 



f(P060) ^ 
f^2020) 



K* 




1()(1500) 



1500) 



^(980) 



f0(98Oi 



f,(400- 
1200 







aj(1640) 



fil520) 



^(1426) Ki(1400) fi 

h(1380) ■ ■ ■ 

fil282) Ki 

^= iq(i270) ■ ■ ■ 

^1 1^(1230) 

hj(1170) 



1(2340) 
1(2300) ■ ■ ■ 

fpi50) 

f^2010) 
f^l950) 

f^lSlO) 
a^l750) /f2 



% ^2 

X(1600) 

f<1565) 

1^1525) 

K*(1430) 

32(1320) 
f2(1270) 



K2 



PDG2000 

TDA 

RPA 



FIG. 5. Scalar, pseudovector and tensor (L = 1) meson spectra. Above 1 GeV, the TDA (dots) 
and RPA (dashes) are essentially identical. 
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FIG. 7. TDA and RPA ground state pion wavefunctions. 
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FIG. 8. D mesons (L = and L = 1). Charmed quark mass is 1200 MeV, light quark is 5 MeV. 
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FIG. 9. Dg mesons (L = and L = 1). Charmed (strange) quark mass is 1200 (150) MeV. 
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FIG. 10. Charmonium system (L = and L = 1). Charmed quark mass is 1200 MeV. 
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